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ABSTRACT 



We present time- resolved optical spectroscopy of V458 Vulpeculae (Nova Vul 2007 
No. 1) spread over a period of 15 months starting 301 days after its discovery. Our data 
reveal radial velocity variations in the He n A5412 and He n A4686 emission lines. A pe- 
riod analysis of the radial velocity curves resulted in a period of 98. 09647±0. 00025 min 
(0.06812255±0.00000017 d) which we identify with the orbital period of the binary sys- 
tem. V458 Vul is therefore the planetary nebula central binary star with the shortest 
period known. We explore the possibility of the system being composed of a relatively 
massive white dwarf (Mi > 1.0 M Q ) accreting matter from a post-asymptotic giant 
k>( ; branch star which produced the planetary nebula observed. In this scenario, the cen- 

tral binary system therefore underwent two common-envelope episodes. A combination 
■ of previous photoionisation modelling of the nebular spectra, post-asymptotic giant 

branch evolutionary tracks and the orbital period favour a mass of M2 ~ 0.6 Mq for 
the donor star. Therefore, the total mass of the system may exceed the Chandrasekhar 
mass, which makes V458 Vul a Type la supernova progenitor candidate. 

Key words: accretion, accretion discs - binaries: close - stars: individual: V458 Vul 
- novae, cataclysmic variables 
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1 INTRODUCTION 

V458 Vul (Nova Vul 2007 No. 1) was discovere d at 9.5 mag- 
nitude on 2007 August 8 (|Nakano et all 120071 ). shortly be- 
fore peaking at V = 8.1. It is classified as a fast nova on 
the basis of its rapid three-magnitude brightness fall from 
maximum within 21 days, indicative of a rela tively massive 
1 Mq) white dwarf. In our first paper IjWesson et al.l 
120081 . hereafter W08) we reported the discovery of a wasp- 



waisted planetary nebula surrounding the r 1 = 18.34 nova 
progenitor, and speculated about the possibility of the cen- 
tral binary star in V458 Vul being composed of a white dwarf 
and a post-asymptotic giant branch (post-AGB) star which 
formed the planetary nebula. However, the lack of an ac- 
curate or bital period pre v ented any further discussion. In 
addition, iGoranskii et alj (|2008h had suggested a tentative 
orbital period of 0.59 d from photometric light curves. In 
an attempt to measure a precise orbital period we started 
a time-resolved spectroscopy campaign searching for the or- 
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Table 1. Log of the observations. 



Date 


Telescope / 


Grating 


Slit width 


Wavelength 


Exp. time 


Time coverage 


Dispersion 


Resolution 




Instrument 




(arcsec) 


range 


(s) 


(h) 


(A pix- 1 ) 


(A) 


2008 Jun 04 


INT/IDS 


R300V 


1.2 


AA3320-8400 


300 


5.7 


1.9 


5.0 


2008 Jun 25 


INT/IDS 


R632V 


1.5 


AA4405-7150 


300 


3.3 


0.9 


2.8 


2008 Jul 04 


WHT/ISIS 


R1200B 


1.0 


AA4929-5644 


300 


7.4 


0.2 


0.7 


2008 Oct 13 


INT/IDS 


R300V 


1.0 


AA3288-8870 


300 


3.6 


1.9 


4.4 


2008 Oct 14 


INT/IDS 


R300V 


1.0 


AA3294-8400 


300 


3.6 


1.9 


1.4 


2008 Nov 11 


WHT/ISIS 


R600B 


1.0 


AA3584-5117 


300 


3.7 


0.9 


1.5 


2009 May 25 


INT/IDS 


R632V 


1.2 


AA4500-6830 


600 


3.2 


0.9 


2.4 


2009 Jul 21 


WHT/ISIS 


R600B 


1.0 


AA4500-4930 


30 


4.4 


0.4 


1.5 




(QUCAM3) 
















2009 Aug 31 


INT/IDS 


R632V 


1.2 


AA4500-6830 


600, 800 


1.0 


0.9 


2.4 



Notes on instrumentation: INT/IDS: 2.5-m Isaac Newton Telescope at Roque de los Muchachos Observatory (ORM), using 
the Intermediate Dispersion Spectrograph with a 2k X 4k pixel E2V CCD; WHT/ISIS: 4.2-m William Herschel Telescope 
(WHT) at ORM, using the Intermediate dispersion Spectrograph and Imaging System with its 2k X 4k pixel E2V CCD. 



bital signature in the radial velocities of the emission lines. 
The results of this campaign are presented in this letter. 



2 OBSERVATIONS AND DATA REDUCTION 

The spectroscopic data were obtained with the Intermedi- 
ate Dispersion Spectrograph (IDS) on the 2.5-m Isaac New- 
ton Telescope (INT) and the Intermediate dispersion Spec- 
trograph and Imaging System (ISIS) on the 4.2-m William 
Herschel Telescope (WHT), both on La Palma. The log of 
spectroscopic observations can be found in Table [T] 

The spectra were reduced using the standard IRAF long- 
slit packages. The one-dimensional spectra were th en ex- 
tracte d using the optimal extraction algorithm of Home 
( 1986). Wavelength calibration was performed in MOLL'fl by 
means of arc lamp spectra frequently taken to guarantee an 
accurate wavelength solution. The spectra were then flux 
calibrated and de-reddened using E(B - V) = 0.63 (W08) 
using MOLLY. For the fast spectroscopy QUCAM3 data we 
averaged in blocks of 10 spectra in order to achieve a proper 
signal-to-noise ratio for radial velocity measurement. 
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Figure 1. Average of 15 spectra taken with INT/IDS on 2008 
June 4 (301 days after the nova explosion). The arrow points 
to the Hell A5412 emission line, which showed significant radial 
velocity shifts while the line emission was still dominated by the 
nova shell. 

binary system at the core of the planetary nebula. This find- 
ing prompted further time-resolved spectroscopy (see Ta- 
ble [TJ in an attempt to measure its orbital period. 



3 EARLY RADIAL VELOCITY VARIABILITY 

The average optical spectrum of V458 Vul taken on 2008 
June 4 (day 301 after the nova explosion) is shown in Fig.[T] 
It is mainly dominated by emission lines of [Nev], [Fevn], 
Hell, and the hydrogen Balmer series. The round-topped 
profile of the Hen A5412 line attracted our attention during 
a first visual inspection of the line shapes. After normal- 
ising the adjacent continuum, we cross-correlated the indi- 
vidual Hen A5412 profiles with a Gaussian template with 
a full- width at half-maximum (FWHM) of 400 km s" 1 . We 
found the radial velocity to vary between ~ —200 and ~ 200 
km s _1 . This radial velocity variation indicated that at least 
one of the components of the Hell A5412 emission forms in a 



1 http://www.warwick.ac.uk/go/trmarsh/software 



4 THE ORBITAL PERIOD OF V458 VUL 

4.1 Period analysis of the radial velocity curves 

We obtained radial velocity curves of the Hen A5412 emis- 
sion line on 2008 June 25, July 4, October 13, and October 
14. By that time no Hen A4686 radial velocity variation 
was detected. This wasn't unexpected since the He n A4686 
line has to thin out (i.e. shed some nova ejecta emission) 
before it starts to present the same phenomenon as the 
weak, optically-thin transition of the Hen A5412 line. The 
Hen A4686 emission showed a clear modulation by Novem- 
ber 2008. Therefore, we also measured radial velocities of 
this much brighter line on 2008 November 11, and 2009 
May 25, July 21, and August 31. Before measuring the ve- 
locities, the spectra were first re-binned to constant veloc- 
ity increments and continuum-normalised. Radial velocities 
were then measured by cross-correlation with a single Gaus- 
sian template. The FWHM of the template used for a given 
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Figure 2. Hen A5412 radial velocity curve obtained with 
WHT/ISIS on 2008 July 4. The velocities were measured by cross- 
correlating the individual profiles with a Gaussian template of 
FWHM = 400 km s _1 . A sine fit to the data results in a period 
of 0.06731 ± 0.00038 d and an amplitude of 115 ± 5 km s _1 . 
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night was adjusted so that the cleanest radial velocity curve 
was obtained, but it always varied between 400 and 1200 
km s _1 . The radial velocity curve of V458 Vul exhibits a 
quasi-sinusoidal modulation. The longest observation (over 
7h, 2008 July 4) covers over five cycles, and a sine fit to 
these data results in a period of 0.06731 ± 0.00038 d and an 
amplitude of 115 ± 5 km s" 1 (Fig. 0. 

In order to refine the orbital period of V458 Vul we sub- 
jected the radial velocity measure ments to a period anal- 
ysis using Schwarzenberg-Czerny's (|Schwarzenberg-Czernv| 
1996) variation of the analysis-of- variance method imple- 
mented as 0RT in MIDAS, which fits periodic orthogonal poly- 
nomials to the phase-folded data. The periodogram calcu- 
lated from the 2008 July 4 data (Fig. O top panel) exhibits a 
strong peak at 14.86 d _1 , consistent with the result from the 
sine fit mentioned above. Next, we analysed the Hen A5412 
radial velocities, which represent about 2/3 of all our ra- 
dial velocity data and were obtained with relatively frequent 
sampling between June and October 2008. The resulting pe- 
riodogram (Fig. [3] middle panel) contains the strongest peak 
at 14.68 d _1 , and the observed alias pattern is consistent 
with the window function resulting from our temporal sam- 
pling. Finally, we analysed the combined Hen A5412 and 
Hell A4686 data, which extend the total baseline spanned 
by our observations to 430 d. The resulting periodogram 
is characterised by a narrow spike at 14.68 d _1 , consistent 
with period determinations of the smalle r radial velocity 
subsets. No signal at the period claimed bv lGoranskii et all 
was found. A sine fit to the whole data set results in 
P = 0.06812255 ± 0.00000017 d or 98.09647 ± 0.00025 min. 
The He n A5412 and He n A4686 velocities folded on the or- 
bital period are shown in Fig. [3] (bottom panel). Our results 
show that the period is coherent for 6341 cycles, suggesting 
it is a fixed clock in the system. We therefore identify this 
period with the orbital period of the binary progenitor of 
nova V458 Vul, which makes it the central binary system of 
a planetary neb ula with the shortest orbital period (see e.g. 
Ide Marccll2009l . for a list). 
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Figure 3. Top three panels: OKI periodograms of the single 
longest observation (see Fig. |2), of all the Hell A5412 radial 
velocities obtained between June and October 2008, and of all 
the Hell A5412 and Hell A4686 radial velocities obtained up 
to September 2009. Bottom panel: Mean-subtracted and phase- 
folded Hell A5412 (green) and Hell A4686 (black) radial veloci- 
ties. The time of zero phase (blue to red crossing) is To (HJD) = 
2454652.52694 ± 8 x 10" 5 . 



4.2 Trailed spectra diagrams 

The long-term evolution of the He n A5412 and He n A4686 
emission lines is shown in Fig. [4] Hen A5412 started to 
reveal the orbital motion of V458 Vul much earlier than 
Hen A4686. By May 2009, Hen A4686 displayed an appar- 
ent orbital signal in the form of a clear S-wave. Note that 
orbital phases were computed relative to the blue-to-red ve- 
locity crossing of this S-wave, which would correspond to 
the standard definition of the orbital phase if the S-wave 
originates on the donor star. The trailed spectra diagram 
of this line also shows high velocity wings extending up to 
~ ±1000 km s _1 . This might indicate the presence of an- 
other emission component apart from the dominant S-wave. 

A deeper look at the May 2009 spectra revealed nar- 
row emission components bluewards of Hen A4686 (see 
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Figure 4. Evolution of the Hen emission lines. Top panel: Hell A5412 trailed spectra diagrams. Bottom panel: Hell A4686 trailed spectra 
diagrams. Black represents emission. No phase binning has been applied, with the exception of the QUCAM data taken on 2009 July 
21. In this last diagram a full cycle has been repeated for clarity. 



Fig. [5]). The first two, counting from Hen A4686, lie at 
rest wavelengths of ~ 4640.6 and ~ 4634.2 A and have 
FWHM ~ 230 km s . These narrow lines, reminiscent of 
the radiation-driven Bowen fluorescence lines used to probe 
the motion of the irradiat ed donor star in X-ray binaries (e.g. 
ISteeghs fc Casaresil2002r ), are in phase with the He n A4686 
S-wave and their radial velocity amplitudes are compara- 
ble within what are necessarily substantial error bars. This 
lends further support to place these S-waves on the irra- 
diated donor star. In V458 Vul, the white dwarf producing 
the nova explosion can provide the EUV radiation needed to 
trigger the process. In fact, two N ill transitions take place 
at 4640.64 and 4634.13 A, very close to the observed lines. 
The other two emissions are likely the Nv doublet lines at 
4603.74 and 4619.97 A. If all these narrow lines originate on 
the irradiated donor star, our adopted phase convention is 
the correct one. 



5 DISCUSSION AND CONCLUSIONS 

In W08 we presented photoionisation modelling of the neb- 
ular spectra obtained before the nova explosion ionised the 
planetary nebula. This implied an ionising source with effec- 
tive temperature T e g ~ 90000 K, luminosity Lboi — 3000 L© 
and radius R ~ 0.23 Rq. In the same paper we showed 
th at, based on the hydroge n-burning evolutionary tracks 
of I Vassiliadis fc Wood! (1994), this requires a core mass of 
0.58 Mq and an age since leaving the AGB consistent with 
our estimated nebular expansion age of 14000 years. 
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Figure 5. Top: Hell A4686 trailed spectra diagram of the May 
2009 run after averaging into 40 phase bins. Contrast has been ad- 
justed to highlight the main S-wave (left) and the narrow compo- 
nents bluewards of Hell A4686 (right). Black represents emission 
and a full cycle has been repeated for clarity. Bottom: Doppler- 
corrected average spectrum of the May 2009 run. 
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The question now is: which of the two stars is the 
progenitor of the planetary nebula? The short orbital pe- 
riod of V458 Vul and the age of its planetary nebula may 
seem at odds if one assumes the system is actually a cat- 
aclysmic variable (CV) which evolved from a much longer 
orbital period b y losing angular momentum due to magnetic 
wind braking (|Verbunt fe Zwaanl Il98ll ; iRappaport et all 
1983T) and radiation of gravitat ional energv jFaulknei 197ll ; 



Paczvhski fe SienkiewicJll98ll ). The time it takes a CV to 
evolve down to an orbital peri od of 98.1 min is of the order of 
a Gyr (IRappaport et alii 19831 ). However, it is possible to get 
a short-period, normal CV within a common envelope (e.g. 
the case of the young pre-CV S PSS J005245.11-005337.2 in 
iRebassa-Mansergas et al1l2008l ) , but producing a nova event 
in this scenario within 14000 years is very unlikely. 

This strengthens the possibility, as suggested by W08, 
of the donor star in V458 Vul being actually an evolved 
star, i.e. a post-AGB star. In such a case, the planetary 
nebula of V458 Vul may have been ejected by the donor star 
instead of the accreting white dwarf after a second common- 
envelope phase. As mentioned above, the post-AGB donor 
would therefore have a mass of 0.58 Mq. 

Theoretical n ova models (e.g. IPrialnik fc Kovet3ll995l ; 
lYaron et al.ll2005l ) agree that a minimum white dwarf mass 
Mi ~ 1 Mq is required to trigger the thermonuclear run- 
away in fast novae like V458 Vul. Observations, although 
scarce, point to a similar value l|Ritter fc Kolbll2003l ) . Hence, 
the total mass of V458 Vul may well be > 1.6 Mq, above 
the critical Chandrasekhar mass, indicating that it may be- 
come a Type la supernova if the white dwarf manages to 
accumulate mass in the presence of nova eruptions. 

Several other systems have been claimed as Type la su- 
pernova progenitors. The subdwarf-B+white dwarf binary 
KPD 1930+2752 is among the best candidat es, but its to- 
tal mass is very close to the critical mass (|Maxted et al.l 
l200d ; lErgmaet all l200ll ; iGeier et al.ll2007l ). The first He 
nova, V445 Puppis, may contain a binary system composed 
of a ma ssive white dwarf ac creting from a helium star com- 
panion jWoudt et al.ll2009l ). The 3.9-h central binary star 
of planetary nebula P NG135.9+55.9 (SBS 115 0+599A) has 
also been put forward (|Tovmassian et al.ll2010l ). In this case, 
a post-AGB star and, presumably, a compact companion 
also amount to a mass just close to the Chandrasekhar limit. 

An obvious objection to our scenario is the fact that the 
post-AGB donor star would have to fill its Roche lobe in or- 
der to sustain mass transfer while it is still contracting. A 
star filling its Roche lobe must obey an orbital peri od-mean 
densit y law, so we used the evolutionary tracks of Blocker 
1 19951 ) in an attempt to find stellar parameters which fit 
both the 98-min orbit of V458 Vul and the results of our 
photoionisation model. We find that a star with an initial 
and final mass of 3 Mq and 0.625 Mq, respectively, on a 
helium burning track at 14000 yr, provides almost perfectly 
the measured effective temperature, luminosity and radius of 
the ionising source. However, steady mass tranfer (i.e. con- 
tact with the Roche lobe) depends on the timescales over 
which orbital angular momentum is lost during the second 
common-envelope phase and contraction of the post-AGB 
donor star take place. Both processes have very short and 
similar timescales (of a few thousands years), making the sit- 
uation very difficult to quantify. Only further spectroscopic 



search for spectral lines from both components of the binary 
system may shed more light onto its dynamics and nature. 

In conclusion, we have solidly measured an orbital pe- 
riod of 98.09647 ± 0.00025 min for V458 Vul. A plausible 
scenario explaining V458 Vul is that of a double common- 
envelope binary system composed of a Mi > 1 Mq white 
dwarf (the accretor) and a M2 ~ 0.6 Mq, post-AGB star 
(the donor) which expelled the planetary nebula 14000 yr 
ago. The total mass of the system may therefore well exceed 
the Chandrasekhar mass which, in addition to its close orbit, 
makes V458 Vul a Type la supernova progenitor candidate. 
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